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ABSTRACT 

It is shown that the diurnal density and temperature variation within the 
thermosphere is generated not only by solar EUV heat input but also by a tidal 
wave from the lower atmosphere which penetrates into the thermosphere and 
predominates the diurnal variations below 250 km altitude, Assuming a tidal 
wave from below which is independent of solar activity and assuming tolar 
EUV heat input to be proportional to the solar activity index F 10 7 , the ob- 
served diurnal density variation and its dependence on solar activity can be ex- 
plained and reproduced quantitatively with help of a two dimensional model be- 
tween 100 and 400 km height. The calculated density agrees with the Jacchia 
model above 300 km and with observations made by Priester, May, King-IIele, 
Taeusch and by Marov below 250 km. 




i. mmmvcwm 

Hliice the launch of the first satellite In 10.17. satellite drag measurements 
and direct in situ observations have rapidly increased our knowledge of the 
neutral density and its temporal variations within the thermosphere. We can 
now distinguish basically between the following regular and sporadic density 
disturbances within the thermosphere (see e.g. the review article by Priester 
et aL, 1007). 


Regular effects; 

la. Diurnal variation (Period i 1 day) 

lb. Solar rotation variation (Period; 27 days) 

la. Semiannual variation (Period; 0.0 years) 

Id. Solar cycle variation (Period; 11 years) 

Sporadic effects; 

2a. Large scale neutral air waves (Quasiperiod; *** 1 hour) 

2b. Geomagnetic activity effect (Quasiperiod; ^ 1 day) 

While the regular disturbances labelled (1) are predictable in occurrence 
and period* the sporadic effects labelled (2) are not. Nevertheless the last two 
effects have simple basic features; a train of harmonic waves in case (2a) and 
an impulse form in case (2b). Thus, they both are accessible to simple model 
calculations , 
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Apart from the disturbances mentioned above there exist of course other 
irregular disturbances without such simple features which we shall not discuss 
in this paper* 

There is general agreement that the effects numbered {la), (lb) and (id) are 
clearly related to solar EUV radiation which is absorbed in the thermosphere 
and transferred into heat of the neutral air* Effect numbered (2b) is connected 
with geomagnetic and ionospheric disturbances within the thermosphere and l p 
thought to be caused by solar corpuscular radiation. Effect numbered (1c) is 
also correlated with geomagnetic activity, but in this case the causal relationship 
is not known, Effect numbered (2a) which is related to large scale travelling 
ionospheric disturbances is caused in or below the base of the thermosphere and 
Is interpreted as free internal gravity wave propagation into the thermosphere 
(Newton at al», 1909). Wo shall exclude this effect entirely In our subsequent 
discussions. 

The effect, best understood, is the diurnal variation (numbered la). The 
theoretical model of Harris and Prlostcr (1902) connected diurnal EUV heat input 
with density and temperature variations within the thermosphere. The difficulties 
of obtaining quantitative agreement between calculated and observed density 
data in the Harrls-Priestor model have been settled recently by Holland ot al. 
(1968). 

As a working hypothesis Jacchia and other people (see e.g, Pricstor et al., 
1967) invented empirical formulae which relate the exospheric temperature of 
the static Jacchia model (Jacchia, 1904) or of the Harris -Pricstor model (C1RA, 
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1 m the vartmiH activity parameters. This method works well between 200 
amt mt km altitude* It fulls however to reproduce some observations below 
200 km e.g, the semlumuutl effect and the diurnal variations. Moreover, this 
method does not provide any Insight Into the physical processes which cause the 
various ci frets* 

Home of the problems more or less unsolved until now are the following! 

1 At what height is the predominant energy deposited for the various 
effects and what form and amount of energy is it? 

2. What are the reasons for the delay between the time of observed density 
maximum for the different effects and the assumed time of maximum 
heating? 

Wo shall try in this papor to give a consistent answer to these questions. 
Part X deals with diurnal and solar cycle variations. In Part IX of this paper we 
shall treat solar rotation variations, semiannual variations and the geomagnetic 
activity effect, 

2. BASIC THEORY 

It has been shown that a two dimensional equatorial model of the thermo- 
sphere can explain fairly well basic data of the diurnal tido within the thermo- 
sphere between 100 and dOO km altitude (Volland and Mayr, 1908a, Volland et al., 
1908), The diurnal variation of density and temperature is generated by solar 
EUV heat input within the thermosphere and by a tidal gravity wave from the 
lower atmosphere which penetrates into the thermosphere (see Figure 1). Heat 
Input q, density and temperature T therefore can be separated in the following 
manner; 
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are amplitudes and times *rt maximum of the disturbance within the thermosphere 
generated by the tidal wave from the lower atmosphere. These terms, too, are 
height dependent and slowly varying with time w ith periods large compared with 
one day. 

u Is the angular frequency of the earth's rotation, 

it* -jp is the local time, 

t is the universal time, y the horizontal distance along the equator, / the 
height and R the earth's radius. 

Ae filUfn and are amplitudes and times of maximum of the total 

diurnal disturbance observed w ithin the thermosphere. 

The calculations are based on perturbation theory which is a sufficient ap- 
proximation below 400 km altitude. Since perturbation theory is linear, the 
variations generated by KUV and by the tidal wave from below within the ther- 
mosphere are completely decoupled from each other. Moreover, the physical 
parameters density and temperature are directly proportional to their respective 
exciting sources Aq puv and Aq ljU ^ , To obtain agreement between the theoretical 
results and the observed density, temperature and horizontal wind only the two 
components of the heat source Aq BUV and Aq ude can be adjusted. Those will be 
chosen so that the density data of Jucehia between 300 and 400 km height and the 
density data of Pries ter et at, (1900), May (1903), Marov (I960, 1908). Taeusch 
et ai, (1908) and of King-Hole and Kingston (1907, 1908) below 200 km are re- 
produced. Then it turns out (Volland and Mayr, 1908a) that the calculated hori- 
zontal equatorial wind at 100 km altitude approaches very closely the amplitude 
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ami phase *if the diurnal wlml derived by ’Jut- front the geumagnet'e S,, 
current! 

We consider only the first harmonies of the diurnal variation, because the 
diurnal wave Is predominant throughout the thermosphere above 1.10 km altitude. 
Thus, we can make use of the following relationship that allows a comparison 
between observations and calculations: If *> and t- are the observed maxi- 

l if *1 i • IH 

mum ami minimum values of a certain physical parameter during one day then 

<5 t p il *» 

f - O * 4 i fj fy_ ^ ^ * 4* >lf« 

p at " ij ti 

HI l* ft * fj I fiffl ft 

Is an excellent approximation m long m the diurnal component predominate#* 

The mean thermospheric density, temperature and molecular weight between 
120 and 400 km altitude have bean chosen from dacehia's static diffusion model 
(Jaechia, 1904) and between 100 and 120 km from (JIRA (I960). By this choice 
wo have eliminated the mean heat input q fi * 

8, TIDAL WAVE FROM BELOW 

A diurnal tidal gravity wave is thermally driven by direct isolation absorp- 
tion by water vapor in the troposphere and by ozone In the stratosphere and 
mesosphere (Lindzen, 1987). Its exciting force — solar radiation within the UV 
and visual light spectrum — does not depend on solar activity. However, its 
transformation Into wave energy and the propagation characteristics of the tidal 
wave depend on the meteorological conditions which vary with the seasons, We 
shall show in Part II of this paper that this seasonal variation of the tidal 
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wave muse# the semiannual efteet of Hu* tin rroospherie el^miiiv. In fh,s 
paper where we #l^al imlv with diurnal ami soiur evade variation* hi* con- 
sider the tidal navi* from la*l*m tin eoimtant ami imlejiendent on solar 
activity# 


Him* the lower imundary oi imp model is at Imm km altitude#* ami thus far 


tram the height range of the generation at the wave* we must treat the diurnal 


tidal wave from helaw m an external energv entree (see ligure 1) and consider 


lilt magnitude ami phase as jmrnmeters* In order la adjust this tidal wave la the 
observed density data of Taeuseh et ah and ot King-llele and iiingston 
(11)117* liM*) as well as to the wind value of Kato (WMl) we ehose the relative 
density at 10ft km altitude nf the tidal gravity wave to la* 


y t id# 
' 0 

p t I 'I* 




0. 12 

20‘ MI lend time J 


f »t /„ 1*** fcr* 


(ft) 


The longitudinal wind component oi this tidal wave at the equator at 100 km 


altitude Is then 


v tiiir 


v t i df* 


16 m see 


V at /, 


local t imi* J 


100 km 


( 7 ) 


which is close to the value derived by Kato (1 twit) from the S t| current. 

* 
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It cun he shown that the time averaged vertical energy flux of that tidal 
wave Is of the order of 10 " 2 erg/cm 2 see at 100 km altitude, This is two orders 
of magnitude smaller than Lindzen's (1007) theory predicts. In Llndzen’s cal- 
culations, dissipation of wave energy due to heat conduction and ton drag has 
been noglee* '» Therefore his energy value Is probably exaggerated at 100 km 
altitude# On the other hand, in that height range the diurnal tidal gravity wave 
belongs to the evanescent modes for which the classical definition of effective 
vertical energy flux 


wave 



p COSh 


w 



on which the above mentioned energy data are based is problematic (Volland, 
1009), A w, ; w and Ap, * are amplitudes and times of maximum of vertical wind 
and pressure components, respecter, of the wave. Note that ;*.r an evanescent 
free internal gravity wave within an atmosphere of zero heat conductivity wo 
have 

r*(7 w "» p ) w/2, thus *% QVe o, 

The tidal wave penetrates Into the thormosphoro and predominates the 
horizontal diurnal wind system below 20 0 km height (Volland and Mayr, 1008a), 
Most of its wave energy is dissipated by heat conduction and ion-neutral colli- 
sions in this height Legion which heats the thermosphere in addition to the EUV 
heating# Wo shall see in Section 8 of this paper that the solar cycle effect pro- 
vides the possibility of estimating the amount of heat due to the dissipated energy 
of the tidal wave within the thermosphere. 
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u iSOLAH cuv heat mm' 


Solar EUV radiation is transferred into heat ot the neutral gas primarily by 

phntodissoefntion of 0> and N\» by Urn recombination and by electron cooling* 

Photodissoeiation, the dominant heat mechanism within the lower thermosphere, 

becomes insignificant at about 100 km and below because at these altitudes almost 

ail of the EUV radiation is attenuated* On the other hand, above 400 km ion 

cooling is the predominant heat source* However, its heating rate there is at *> 

least two orders of magnitude smaller than the healing rate at altitudes between 

100 and 150 km* Because of these reasons we confine our thermospheric model 

to a height range between 100 and 400 km and assume that the variable part of the 

solar spectrum — namely EUV radiation - is deposited mainly In this altitude 

range. 

* 

The height and time dependence of the various heating mechanisms is not 
very well known. Harris and Prioster (1902) assumed that the heat Input is 
proportional to the Ion production rate. Their EUV heating rate for CIRA 
model 4 (F 125) can be approximated by the following analytical expression 







^EUV 


■»f *“Z0 ) %uv 


r EUV 12° 0 local time 


(8) 
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Herts ‘ E evv ium erg/(em 2 Hee) is the diurnal component of the total effective 
heat input into the thermosphere above 100 Urn* z t# 100 km is the lower boundary 
of our model and H pyv fd* km is an empirical seale height. 

We shall adopt a heat input of the general form of Equation (H). However, 
because of our lack of knowledge we consider " k K!?v and H Pnv as free parameters 
in order to adjust our calculations to Jacobin's density data between 000 and 
400 km altitude. For simplicity we consider the scale height II Bl , v as a constant 
independent on solar activity. The total diurnal EUV heat input is supposed to 
be directly proportional to the solar activity factor F (F is the 10.7-em flux in 
units of io~ 22 Watt/m 2 Hz smoothed over two or three solar rotations). 


The optimal parameters which for reproducing Jacchla's data in 
range between 000 and 400 km are 


MS 


EUV 


0.27 


F 


125 


c*rg / (cm 2 sc‘c) 


^buv 


( 0 ) 


Wo shall use those parameters in the following calculations# 

Our scale height H £uv is larger by a factor of two than the value derived 
from Harris and Pries tor's (1902) calculations. This is not surprising because 
we shall see in section 6 that the EUV heat input significantly influences the ther- 
mospheric density variations only above 250 km altitudo. The heating of the 
thermosphere at these heights is provided mostly by the plasma of the ionosphere. 
Thus , the height distribution of the heating rate not only depends on the neutral 
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density but al8o on the density and temperature distribution of the electron*. 

Our scale height II Frv approximates the Ki V heat distribution between alsmt 
aoo and ion km altitude. The real heating rate certainly is much more complex 
than Equation (*) exhibits, 

f>, COLLISION NUMMSH 

Collisions between neutrals and ions play a decisive role in the dynamics ji 
the thermosphere at 1*' 2 -layer heights* There, collisions between oxygen-atoms 
and oxygen-ions predominate, Dalgarno (loot) gave the formula for eoiiisions 
between oxygen-atoms and -ions: 

- ur l(i n i see* 1 (10) 

where the factor / varies between 5 and 10 at ion temperatures bevween T, 

400 ‘ and 2000 K, and N, (in enr 3 ) is the number density of O-ions. Near the 
equator the real height of the F 2 -layer Is between 000 and 400 km, We approxi- 
mate the time averaged collision number c by the analytical expression 


1 

3 



\f o 

mux ' 


l> f ? - z wo x ) * 


( 11 ) 


a# 


where max r 4 x 10" 4 sec* 1 is tiie maximum collision number at the maximum 
height z mttX - 200 km for F r; 125, and b - 7 x 10~ s km* 2 determines the 
height dependence of ?*. The bracket on the right side of Equation (11) allows 
for a variation of v by a factor of two between solar minimum and solar maxi- 
mum, According to Equation (10) the value of v in Equation (11) for F 120 
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corresponds to an inn density amt an electron pkrnma irequeney f t , respectively, 
of N\ 5,5 * lo s tm < and f p U.7 MI!/. 


Equation (Id) in oversimplified at K -layer heights, Tin re, the collision cross 
section decreases due to the change in the major constituents ol plasma and neutral 
air. However, because ol the small Ion number density at this height range the in- 
fluence ol ion drug on the diurnal wave becomes small. It can be shown that tin 4 


exact altitude profile of , determined by the exponential factor U In Equation (11), 
is not of great significance for the diurnal dynamics, However, tin* absolute value 
of maximum collision number ri((Jt within tin* maximum F2 -layer has an important 
influence on the time response of the thermospheric density with respect to solar 
KUV heating. This is shown in Figure a where the time of maximum r rrv of 
the diurnal density variation at arm km altitude has been calculated for values 
of i* max varying from lo * n to lft* 2 see 5 , We note immediately that the observed 
time delay of 2 hours in the density variation (maximum at lT m local time) with 
respect to the KUV heat input (maximum at 12 ,m local time) corresponds Just to 
the collision number derived from Equation (10) and from the observed time 


averaged maximum electron density of the equatorial F2 -layer (see upper scale 


in Figure 2 where the ion number density is scaled). 


The observed time lag of two hours of the thermospheric density is there- 
fore the natural time response of the thermosphere to the KUV heat input. The 
dependence of this response time on the collision number results from the strong 
influence of . on the horizontal winds. The maximum longitudinal wind speed 
Av Et?v of the diurnal wave at 350 km altitude is plotted versus in Figure 2 
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ah dashed line. The collisions act like a “resistance” within an oscillator system. 
On the other hand, the horizontal winds determine the convective heat capacity of 
the thermospheric “enndensor”. An increase of the eoilishm number leads to 
a decrease ot the horizontal wind speed, thus to an increase of the convective 


heat capacity and an increase of the time lag. For small horizontal velocities 
the time lag becomes 5 hours (corresponding to a density maximum at I7 m local 
time). This is just Harris and Dniester's (1902) result in their one dimensional 
vertical model in which horizontal winds have been negleeted. A decrease of the 
eollision number gives rise to a phase advance of the density maximum with re- 
spect to the KU\ heat source which becomes 00 (or 0 hours) at small collision 
numbers. This is the direct result of the decreasing resistance and of the de- 
creasing convective heat capacity within the thermospheric "circuit system”. 


The relative amplitude A, mv / n) is not influenced very much by the collision 
number nilH . 

0* DIURNAL AND SOLAR CYCLE VARIATIONS 

Using the assumptions about the energy input made in section 2 to 5 and 
using Jaeehia‘8 static diffusion model for the mean values [Jacchia (1904)] 
our thormosphexTc model of the diurnal variations over the period of one solar 
cycle is uniquely determined. Figure 3 shows the relative amplitude of density 


and temperatiue variations versus height due to the Internal EUV heat input 


[Equations (8) and (9)] for four different solar activity factors F - 02.5, 125, 
187,5 and 250, equivalent to Jacchia 1 s exospheric temperatures of - 750% 
1000°, 1250^ and 1500° K, Though the heat input increases with increasing F, ine 
relative amplitudes decrease at altitudes above 200 km. This results from the 
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fact that mean densitv ( „ ami mean temperature 'i,, Increase wirt* with ¥ 

Hum the diurnal variation. However, tin* absolute amplitudes |( v , ami 1, p v at 
-tint km altitude, plotted vermis ¥ in Figure 1, increase with K. Below 2uo km 
the strong heat input causes a minor increase of the relative amplitude with In- 
creasing solar activity. However, the relative density *arialiou due to FFV is 
smaller than 2 there and In el no significance. 

Figure .1 shows the relative amplitudes oi density and temperature for the 
contributions due to the tidal wave from below , The tidal wave propagates freely 
through the thermosphere* Between loo and ion Urn altitude about fK of its wave 
energy is dissipated by heat eonduetion and ion drag* Therefore the relative 
density amplitude remains roughly constant with height while the relative tem- 
perature amplitude decreases between ir>a and a on km altitude* The decrease of 
the relative amplitudes with increasing solar activity is again the result of the 
strong increase of the* mean values with F. In Figure -1 we notice that the absolute 
density variation generated by the tidal wave increases with F, while the tempera- 


ture variation remains constant. This results from the changing propagation eon 
ditions with in the thermosphere with solar activity because according to our as- 
sumptions the tidal wave input at 100 km does not depend on solar activity. 


The sum of both contributions is plotted in Figure 0 versus height for the 
four exospheric temperature values 1\ . The dashed lines in Figure 0 are 
Jacchia’s equivalent day to night ratios converted into relative amplitudes by 
Equation (3). The full and open circles are N 2 measurements by Tucuseh et al» 
(1908) and satellite drag measurements by King-IIole and Kingston (1907, 1908) 
made at moderate solar activity (F 100 - 100). Wo repeat that the four 
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aval! <h!e iiurnmHerg *4 tlu* |«*t heat uoureeH m mir iiunld have bt en chosen 
such that 4urehiu*s data beuu en unit amt ion kin am! iuettHelt’H ami King-Bele*# 
data below 2»M km at F 12 * could In* n produced. 


‘Hu* results show that in view ot »nu* highly idealized model tlu* agreement 
between Juerhlst'H dennilv data am! our euleuiatioms ih exeellent above Him km* 
Hr* dependence of tiu* relative density amplitude* »m Holur activity In well repro* 
dueed, The largest dliTemiecH appear at highest solar activity. It is of course 
possible to obtain even better agreement by allowing an increase of the scale 
height II H v with solar activity oi the id V heal energy rate* in liquation ts). 

Below 2.10 km altitude we notice a strong deviation oi our density data from 
Jaeehia's densities. This is due to the iaet that Jaeehia forced his density varia- 
tions to become zero at 120 km altitude while In our calculations the tidal wave 
predominates at this height range. At 200 km altitude Priesler et al* (1000) and 
May (1000) obseived relative density variations ol 10 , in 10.1* during higher 
solar activity while King-Hole and Quinn (100.1) and Marov (190.1) during low 
solar activity in 1000 gave day to night ratios ol* i 1,7 to 1,0, These data are 
correctly reproduced in our model, Our temperature amplitudes are systemati- 
cally lower than Jaeehia’s data above 200 km altitude. The dilTorenee increases 
with increasing solar activity. This discrepancy will be explained in the next 
section. 


Figure 7 


gives the time of maximum oi density and temperature variations 


due to the two heat sources versus height lor *i\ 1000 K (F 12,1). We notice 

that the KUV contribution (dashed lines) has the time ol maximum at 14 0f ' local 


time for density and temperature above 200 km 


The tidal density maximum 



v4*i*S( h from 20 ' '*ut inu km to l t'“ M at luu Km altitude while tin* itM tempi rotutv 
in nt iu m |im km giving rim* to a phase difference **t i hour* witii mspeet to 
tin* KI"V generated temperature (dash-dotted lines*, The time ot maximum mi 
tiu* total densitv variation (lull lines) remains within d u time above km km 
which in in agreement with observations oi Murov Ulug«iiele and liingston 

(IUU7) and Tuettsch et ai* Ainu the time omimmv* Mtul temperature 

variation which in at local time above Jf*0 Km occurs at Th rime in reamm* 
able agreement with observations of Tut useh et ah (iMn^n 

The times of maximum are very insensitive to solar activity and show 
deviations Irom the data plotted in Figure T which generally are not larger than 
no minutes* 


In the whole altitude range above 100 km there is no isopyenie layer where 


the density variations become zero, 


Figure n shows the diurnal horizontal and vertical wind speed plotted versus 


height and calculated for different solar activity factors F, The horizontal wind 
(positive in the east direction) decreases with increasing solar activity in the entire 
height range while the vertical wind is not influenced very much by solar activity* 
The horizontal wind has a relative minimum near 000 km altitude, which results 
from the maximum ion drag near the F2 -layer maximum, The KIT and the tidal 


contributions of the winds have been plotted in Figure h as dashed and as dash- 
dotted vines respectively only for F 120 (f^ 1000 K) In order to show their 

typical behavior. We notice that the tidal component of the winds predominates 
below 000 km altitude. The horizontal wind responsible for the S f| current at 
E -layer heights is entirely due to the tidal wave from below. 
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hi Figure U |i*e time Ml it* t siiMiiii* m| the »hu* n nt wind eMinfiMNeiiUt uml m! 
the total wind lieen height tor l 12*. Attaint Hit h time ot 

maximum i*onh insigniticmiih iiiilmmedm solar ttetivitv. 

the squares In Figure# *u mu! fin Indicate Kuto** Oh Wf wind of the diurnal 
eompment of the F, current at tin comitor* ihev ill well with the eaienlaied 
values. 


The horizontal wind strongly depend# on the collision number * tin height 
profile of wan given In Fqunthm n it* Therefore we expect a very scnsitP e 
reaction of the horizontal wind m changes in the •profile. This m in hie t the 
case. For example, the calculation# by Tolland et ah (IhfW) were towed on a 
-profile which increased evpimentiaih with height, The horizontal wind there- 
fore did not show the relative minimum a# 1 igure *a indicates, However, the 
wind speeds are of the game order of magnitude. Ainu the time* of maximum 
are nearly the name. Moreover, all other parameters like vertieai wind, density 
and temperature are not significantly Influenced by that change in the -profile. 
This show s that the horizontal wind essentially exerts only a second order effect 
on the diurnal variations* Its integral amplitude can shift the phase* of the* density 
maximum (see Figure 2} hut does not influence very much the amplitude of the 
density, 

7* LIMITATIONS OF TIIK MODEL 


Wo have seen in Figure a that our temperature amplitudes are systematically 
lower than Jacchia's data above 200 km altitude. The reason for this discrepancy 
is twofold. 



Hr**, agmiuent lielwi en Imth model temperatures can he 4 k s|n-rUnl only it 
f^tiii density profile >* are the same* rhis i* not the ease during high solar 
aetivfly ml at altitudes Mm mi km. Our density scale height in smaller i\m\ 
daeehia*#* density scale height fur F mu Therefore our temperature variutt*m 
mum In* smaller than Juechlu 1 * values In accordance w 1th the results In Figure h* 


tfecond, in our calculations we neglected the time variation of molecular 
weight which In due to changing composition from day to night. The barometric 
height formula in an excellent approximation even for the diurnul variation*. 

Our temperature values should therefore agree with Jaeehia 1 # data it hoth density 
profiles are the same. From Figure ha we notice that both density profiles agree 
at T* 7o0 and l mm alsive mm km altitude. JaeehhFs temperature gradient 
is nearly zero almve 800 km height, Our temperature gradient is also small and 
con he taken as zero. Therefore the barometric height formula becomes 


where 



II 


RT 

Mg 



is the density scale height assumed to be equal in JacchkFs model and in ours. 
The diurnal variation of the scale height is then 



(id) 


* 


where (AT/T) r is the relative temperature amplitude of the Jacchia model and 
(AT/T) V is our temperature amplitude because wo neglected the diurnal variation 


in 


*#t till* lit n weight VI* I he dilleieMr irttween ImII* letiijier tfUi'e amplitudes 

in tiue I** the dimiml change >*? M* At im* km altitude and tor I , I mm K It in 


** III O i*Mt >* »h 7 

which in u*rv near i»* n«r value calculated in Figure ah* We MiiouUI however 
Iii-nr in mind fl* *t tin* absolute error in tlu* temperature amplitude lor l\ 
l mm K in nut larger limn ;m K* 




.urmiimg t»* the J.ii * nia model 


| liei eioi e 




vi 


The maximum dliiereme lietweett dueehiiFx and our temperatures Is DU K 
ut luu km altitude during mavlman solar activity (see Figure I). In Figure 10* 
maximum ami minimum temperature* at iim km altitude are plotted versus solar 
activity factor F together with dacehiu 1 * doy ami night exospheric temperatures 
(dashed lines)* We notice that the (inference Mween daeehla's temperatures 
and ours never exceeds d of the total temperature which is certainly within the 
range oi accuracy of any direct temperature observation* 

It is possible to remove the error In the temperature variation of our model 
which is due to the neglect of the temporal composition changes* This cun be 
done by an iteration process* Furthermore, for the expressions of the KFV 
heating rate i Ktjuution (x)| and of the collision number [Ktjuuthm (13)1 simple 
analytical functions have been chosen deliberately in order to iind out their basic 



inliuetm* on Uu* diurnal dynamics, Improved ami mure complete observational 
data especially in the height range below ;*im Km will certainly modify the numeri- 
eal values of these parameters amt ot the tidal wave input* Our model Is flexible 
enough to allow tor such modifleatlonK, it in In fact capable of determining the 
height distribution of the solar KV V heat sourer and oi tin* collision number if 
appropriate density and wind data are available* In this eonneetlon it in perhaps 
of interest that the model calculations for the entire solar cycle on whieh the 
data plotted in Figures a to in are based upon have been performed in the time 
of one minute during a single computer run on the HIM 00 0/7.3, 

Our two dimensional model of course is not able to explain latitudinal effects. 
The global behavior of the diurnal dynamics must be treated with help of a three 
dimensional spherical model (Volland and Mityr, lOttgb). That the two dimen- 
sional model nevertheless describes rather well the low latitude behavior of the 
diurnal tidal variations of the thermosphere is clue to the fact that the diurnal 
thermospheric tides are predominated by the spherical function P u which is an 
eigenfunction of the spherical problem. P u sin * ( is the co-latitude) varies 

only slightly with at low latitudes and therefore can be treated as a constant 
there, 

A word must bo said about the validity of perturbation theory which wc used 
entirely in our calculations. We state that perturbation theory is an excellent 
approximation as long as the relative amplitudes of the perturbated values do 
not exceed the number 0.5. This condition is fulfilled in our calculations within 
the height range between 100 and 400 km for all parameters with the exception 
of V fUum '* 0 at low solar activity ( T. v * 750 U K; sec Figure 0). The error due 



to perturbation theory can host fie estimated In the equation of state. The relative 
pressure amplitude is 

Ap A, AT AM /, AT /14 . 

n ^ ? = t ==3 * m 1 / Cf‘l) 

Pm ; o 4 ii m a f o ‘o 

Apart from the variation of molecular weight M which we already disem ;ed we 
neglected the last term on the right side of Equation (14) In our perturbation 
theory. At moderate solar activity ( 1000 K) and at 400 km altitude this 
term lias the value (see Figure 0) 

A. ■ A q 1 

= 0.5* 0.1 0,05. 

■ 0 A 0 

t 

However, this error of fi% does not appear in the diurnal component but is trans- 
ferred into the semidiurnal component and into the mean value p 0 • The semi- 
diurnal component likewise transfers part of its energy content into the diurnal 
component, As long as the semidiurnal component is small, our perturbation 
theory is valid within the limits mentioned above, It gives of course only the 
first harmonic or the diurnal component of the diurnal disturbance. Fortunately 
the semidiurnal component is in fact small at least above 100 km altitude as the 
observations confirm, This smallness is caused by the filter action of the 
thermosphere which suppresses higher harmonics of the propagation modes 
(Volland, 19(59), 

8. CONCLUSION 

By treating a two dimensional model of the thermosphere at low latitudes it 
has been shown that the diurnal density disturbance is gonerated by solar EUV 
heat input within the thermosphere and by a tidal wave from the lower atmosphere 
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which penetrates into the thermosphere and predominates the density variations 
below Si£iO km altitude. Assuming a constant tidal wave input and a KTV heat in- 


put proportional to the solar uetivity factor l* lu 7 the observed diurnal density 
variations and their dependence on solar activity can be quantitatively reproduced 


between 100 and 100 km height. The calculated density agrees with the Taeehia 


model (Jaeehia, 1001) above 000 km and with measurements made by Prlester 


ct nU (1900), May (1900), King-llele et al» (1000, 1007. 1 90s), by Tneuseh ct aL 
(100H) and by Marov (1000, 1007) below 2f*0 km. 


The time of maximum of the density variation is 14* m local time above 200 km. 
It Is the natural response time ol‘ the thermosphere due to the KUV heat input, 
Below 200 km, where the tidal wave predominates, the time of maximum shifts 
toward 20 mt at 100 km altitude. 


The calculated temperature amplitudes are smaller than Jaeehia’s data above 
200 km altitude. These differences increase with solar activity. The reason for 
this discrepancy is partly due to the difference in the density profiles of both 
models. However a relative error of about iY< must be attributed to the neglect 
of the diurnal variation of molecular weight in our calculations. The absolute 


difference between Jaeehia's temperatures and ours never exceeds O'* of the 


total exospheric temperature. 


Jucchia's (1004) empirical formula for the exospheric day and night time 
temperatures and ihefr dependence on solar activity (the dashed lines in Figure 10) 
can be interpreted as follows,* A contribution of a tidal wave from below which is 
independent of solar activity is superimposed on a contribution due to direct solar 
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I TV heat Input# The uiiiplltiuk‘ m this Id \ gtmToud wnu* mrmtses with solar 
activity and causes the increase (it the exospheric temperature inriuthm with F, 
while the temperature variation due to the tidal wave remain* constant* The wave 
energy of the tidal wave 18 nearly completely dissipated within the thermosphere* 
This dissipated energy contributes to the time averaged heat input into the ther- 
mosphere* We suggest that it is responsible for the residual exospheric* tempera- 
ture of about 800 K at zero solar activity I* o (see Figure 10). From Jucchla'a 
temperature profiles one can estimate that a heat input of about 0*1 erg/ene'see 
above 100 km altitude is necessary to maintain an exospheric temperature of 


800 K. Therefore , we estimate that the tidal wave at loo km altitude has a wave 
energy of 0.1 erg/em 2 soe. That amount oC heat energy contributes significantly 
to the thermospheric heat balance as can be seen by comparing It with the mean 
KUV heat input which varies from about 0,2 to 1 erg/cm *sec between low and 
high solar activity, 


Our two dimensional model can explain basic features of the diurnal varia- 
tions within the thermosphere. However, no attempt was made to explain latitu- 
dinal variations or composition effects like the Helium bulge, Our model never- 
theless describes rather well the low latitude behavior of the thermosphere 
which is due to the fact that the diurnal thermospheric tides are predominated 
by the spherical function P u which is the eigenfunction of the spherical problem 
(Volland and Mayr, lOOSb). 

In Part II of tills paper wc shall discuss the 27 day variation, the semiannual 
effect and the geomagnetic activity effect. 
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Figure? 1, Block diagram showing the two energy sources which generate diurnal 
density disturbances within the thermosphere. 
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Figure 7. Time of maximum of the diurnal variations of density and temperature 
calculated for a thermospheric model of exospheric temperature of T - 
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Figure 9. Time of maximum <4 the diurnal variations of h oriental arid vertical Hands, 
calculated for a thermospheric model of exospheric temperature of T x ” 1300 K, 



Figure JO* Calculated maximum and minimum temperatures at 400 km altitude versus solar activity 
factor F. The dashed lines give the exospheric temperatures of the Jacchia model at day and at 




